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^ , We compute properties of the interface of the 3-dimensionaI Ising model for a wide range of temperatures and 

O ' for interface extensions up to 64 by 64. The interface tension a is obtained by integrating the surface energy 
I density over the inverse temperature l3. The surface stiffness coefficient k is determined. We also study universal 
quantities like £_^a and The behavior of the interfacial width on lattices up to 512 x 512 x 27 is also 

investigated. 
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Here we present a numerical study of 3-d Ising 
interfacial properties with focus on the precise 
determination of interface tension and interface 
stiffness over a wide range of temperatures. A 
detailed exposition of our work can be found in 



1. THE MODEL 

We consider a simple cubic lattice with exten- 
sion L in X- and y-direction and with extension 
T = 2D -I- 1 in z-direction. The lattice sites 
i = {ix,iy,iz) have integer coordinates, and the 
z-coordinate runs from —D to -\-D. The Ising 
model is defined by the partition function 



Z 



exp(-/3i/), 

{^.=±1} 



where 



H = — ^ ^ kijCJiaj 

<ij> 



(1) 
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The lattice becomes a torus by imposing geomet- 
rical periodic boundary conditions in all three di- 
rections. For the Ising spins Ui we will use two 
different boundary conditions: Periodic bound- 
ary conditions are defined by letting kij = 1 for 
all links <i,j>. To define antiperiodic boundary 



conditions in z-direction, we also set h 



1 with 



the exception of the links connecting the upper- 
most plane {z = +D) with the lowermost plane 
(z = ^D). These links carry an antiferromag- 



netic factor kij = — 1. 



2. INTERFACIAL PROPERTIES 

We adopt the following definition of the in- 
terfacial width: A magnetization profile for lat- 
tice planes perpendicular to the z-direction is de- 
fined by M(zz) = - (Ti. The antiperiodic 
boundary condition allows us to shift the configu- 
ration in z-direction such that the interface comes 
close to iz = 0. We introduce an auxiliary coordi- 
nate z that assumes half-integer values (labelling 
positions between adjacent lattice layers perpen- 
dicular to the z-direction). Following a nor- 
malized magnetization gradient is defined as 



M{z+^)-M{z-^) 
M{D) - M{-D) 



(3) 



For a given spin configuration, the position of the 
interface is defined as ^ p{z)z. The square of the 
interface width is then defined lH,B| as 



w^^{Y,p{z)z^-{Y.p{z)zy). 



(4) 
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Especially on small lattices, fluctuations in the 
two bulk phases can deteriorate the results. Due 
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to bubbles, p{z) can be accidentially large even 
far away from the interface position. We therefore 
also consider the interface width measured on the 
configurations after a removal of all bubbles. 

We now turn to the definition of the interface 
tension. Let us assume that there is exactly one 
interface in the system with antiperiodic bound- 
ary conditions, and that there are no interfaces 
in the periodic system. Then the interfacial free 
energy is Fg = Fa.p. — Fp, + In T, where Fa.p. (Fp.) 
is the free energy {— — In Z) of the system with 
antiperiodic (periodic) boundary in z-direction. 
For a discussion of the case of several interfaces 
see [0. The interface tension is defined as the 
limit a = limi^tx) Fajl? . 

In general one has no direct access to the parti- 
tion function in Monte Carlo simulations, except 
for not too large systems, cf. [Q. Note, however, 
that 



2-dimensional Gaussian model with Hamiltonian 



d(3 



{H)a.p. — {H)p. = Eg 



By integration one obtains 
FM=FAf3o)+ / dp'Egip'), 



(5) 



(6) 



where /3o is arbitrary. Our approach is to compute 
by Monte Carlo simulation the surface energy for 
/3-values ranging from the critical region around 
(3c (= 0.221652(3) g) up to /3o = 0.6. For large 
(3 we can employ a low temperature expansion by 
Weeks et al. (published in an article by Shaw and 
Fisher |^) to obtain Fs{(3q). Note that it is also 
possible to start the integration at small (3 
(To the best of our knowledge, the integration 
method to obtain surface free energies was first 
used by Biirkner and Stauffer Q]). 

In the theory of rough interfaces the surface 
stiffness coefficient k plays an important role. If 
by suitably chosen boundary conditions the low 
temperature interface is forced to make an angle 
9 with e.g. the x-axis, one defines (see e.g. [pi). 



(7) 



Capillary wave theory says that the long distance 
properties of the interface should by encoded in a 



Ho 



2/3eff 



(8) 



and /3cff = Long distance properties are 

most systematically studied via the block spin 
renormalization group. For the Gaussian model 
defined through eq. (^) one defines block spins 
$/ as averages over cubic blocks / of size Lg, 
i.e., $/ = Lg^ J2iei define the quantities 



a: 



(0) 



(9) 



where / and J are nearest neighbors in the block 
lattice for i = 1, and next-to-nearest neighbors for 
i = 2. I is the extension of the block lattice, i.e. 
I = L/Lb- For the Gaussian model, the ^'s can 
be computed exactly. For the Ising model, block 
spin "height variables" hj are defined as follows: 
The blocks / are sets that are quadratic in x ^ y 
direction with extension Lb ^ Lb and that ex- 
tend through the whole lattice in z-dircction. For 
every block, a magnetization profile and an inter- 
face position can be determined exactly as in the 
case of the full lattice. We define hj = interface 
position in block /. Note that the blocked height 
variables can also be defined "with and without 
bubbles" . Blocked observables for the Ising inter- 
face are introduced analogously to eq. (^ , and are 
denoted by For a rough Ising interface, 

we define an effective coupling PcS as follows: 



/3cff = lini 



A 



(Ising) 



(0) 



(10) 



Of course, we expect that the so defined /?off does 
not depend on i or L 

3. MONTE CARLO RESULTS 

We did simulations with antiperiodic bound- 
ary conditions in z-direction on lattices with L = 
8,16,32,64 for /3-values ranging from the bulk 
critical region up to /3 = 0.6. For many /3-values 
we made runs with different D to control the ef- 
fects of a finite thickness of the lattice. In total, 
we made more than 250 different simulations with 
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antiperiodic boundary conditions. Typically, we 
made 10000 measurements of several quantities, 
separated always by 8 cluster updates with the 
Hasenbusch-Meyer cluster algorithm for Ising in- 
terfaces The simulations supplied us with a 
sufficiently dense grid of /3-values for the energies 

For most of the values, we fortunately did not 
have to do extra simulations to access the ener- 
gies with periodic boundary Ep, . Instead we used 
the diagonal Fade approximation of the low tem- 
perature series by Bhanot et al. |10|. For smaller 
/3-values we used the cluster Monte Carlo method 
to determine Ep, . 

In order to do the integration over /3 we interpo- 
lated the data with cubic splines. The integration 
was then started at /3 = 0.6, where the integra- 
tion constant can be safely taken from the low 
temperature series. 

Already for moderate surface extension L, the 
surface free energy was found to behave with very 
good precision like Fg = Cs + a' L^. It was there- 
fore natural to identify the coefficient a' with the 
surface tension a. 

The results for the free energies for L — 8, 
16, 32 and 64 were then used to make fits with 
Fs = Cs + (J L'^ in order to obtain estimates for 
the surface tension a. In table || we display a 
few of our results. Our results show a significant 
deviation from a prediction by Shaw and Fisher 
|6| based on an analysis of the low temperature 
series. A detailed comparison will be published 
elsewhere [Q. 

We fitted our results for a to the critical law 
a = U{)t^ , using both of the two definitions 
t = 1 - Pc/ P and t = (3/ (3c — 1- We also var- 
ied the interval over which the /3-dependence of a 
was fitted. The fits were always done using four 
different /3-values. The results based on the two 
different definitions of t were statistically incom- 
patible, showing that one is still not close enough 
to criticality. However, taking systematic effects 
into account, we consider our results consistent 
with fi Ri 1.26. The results for the critical am- 
plitude (To show even stronger dependency on the 
type of the fit, and we can not say very much 
more than that Indo is probably something be- 
tween 0.2 and 0.4. 



In order to study the behavior of the product 
^^cr we tried to determine the correlation length ^ 
from the simulations of the system with periodic 
boundary conditions. We defined ^ via the de- 
cay of the connected 2-point correlation function 
of the absolute value of the time-slice magnetiza- 
tion. Our results are nicely consistent with low 
temperature series [ pT| . 

We measured the block spin correlation func- 



tions Afi"^'^' and studied the quantities /3^g(A) = 

Aff"'^\LB = \)/A[JiLB = oo) (measured on 
the configurations with the bubbles removed). 
The values for the different i, I and A's turned 
out to be fairly consistent within the statistical 
accuracy. Our estimates for (3cS were then deter- 
mined by averaging over /3^g with i = 1 and i — 2. 
Some of our results are shown in table |^. 

In figure |^, we show our results for two com- 
bined quantities, namely ^^ct and ^^k. In the 
product ^^cr, the exponents n and —2v of the re- 
duced temperature t should cancel, and we expect 
that this product should be fairly constant in a 
neighborhood of the critical point. The full line 
in the figure was obtained by combining our a's 
from the integration method with the correlation 
lengths as obtained from the Fade. 

Since we do not know the error of the Fade ap- 
proximation of the low temperature series we base 
our error estimate for this quantity on our error 
bars for the measured correlation length and on 
the statistical errors on the surface tension a. We 
estimate the relative precision of our results for 
^^(7 to be around 5 per cent for the smaller /3's, 
certainly better in the large /3 region. This takes 
into account statistical errors only. There might 
also be systematic errors (due to too small L's) 
in the surface tension close to the critical point. 
They might be responsible for another 5 percent 
relative uncertainty. The points with error bars in 
the figure show the product ^^k. The plot shows 
that in the critical limit the surface stiffness be- 
comes the same as the interface tension. This 
is a consequence of the restoration of rotational 
symmetry at the bulk critical point. Using our 
results for both ^■^ct and ^^k, we estimate that in 
the limit /3 — > /3c both quantities have the limiting 
value R- = 0.90(5). 
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Table 1 



A few results for a and for /3off 



/9 


0.2275 


0.2327 


0.2391 


0.3000 


0.40236 


a 


0.0146(1) 


0.0319(1) 


0.0555(1) 


0.30284(8) 


0.67988(6) 




0.240 


0.275 


0.330 


0.35 


0.37 


Pes 


16.7(5) 


4.65(4) 


1.93(2) 


1.52(2) 


1.20(2) 



Figure 1. Results for f^cr and 



In the theory of critical wetting, the quantity 
— l/(47rK^^) plays an important role. In 
a paper of Fisher and Wen |l^ this quantity is 
estimated over a wide range of temperatures. A 
comparison of our data with their theoretical pre- 
diction will be published elsewhere |p^ . 

When (3 approaches the roughening coupling 
Pr, Kosterlitz-Thouless theory states that f3cS ^ 
2/tt. Using the estimates Pr = 0.4074(3) (m), 
and ^(/3cff) — 0.3163 (from the Fade that here 
certainly is reliable), we find a "KT value" of 
which is 0.1572. 

In order to demonstrate the efficiency of the 
Hasenbusch-Meyer algorithm we redid the sur- 
face width computation of Mon et al. |Q at /3 = 
/3c/0.8 = 0.2771 on lattices of size L x L x 27, 
with L = 32, 64, ... , 512. We performed fits of 
the data using the ansatz 



that is motivated by Kosterlitz-Thouless theory of 
a rough interface. When the L ~ 32 data were ex- 
cluded, the fits where quite convincing. From our 
analysis we arrive at an estimate /3cS — 4.3(2). 
This result is nicely consistent with the (3cS as 
obtained from the renormalization group quanti- 
ties A. 

M.H. would like to thank the Deutsche 
Forschungsgemeinschaft for financial support. 
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